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SUMMARY 


The  purpose  of  this  research  was  to  develop  an  effective  program  for  comput- 
ing the  resistance  to  ground  of  simple  electrodes  carrying  direct  current.  A model 
of  horizontal  earth  layers  and  parameters  applicable  to  arctic  engineering  was  em- 
ployed. The  report  begins  with  a brief  introduction  to  the  theory  of  d-c  earth 
conduction  and  presents  the  definition  of  resistance  to  ground.  Tire  development 
of  the  specific  electrode  cases  of  a horizontal  straight  wire  and  a vertically  driven 
rod  are  then  discussed.  The  explicit  developments  of  the  mathematical  solutions 
are  not  presented  because  they  may  be  found  in  the  literature.  However,  the 
solutions  themselves  are  presented  and  they  arc  then  numerically  integrated  us- 
ing earth  parameters  based  on  previous  studies  of  permafrost  resistivity.  The  re- 
sults of  the  simpler  cases  involving  homogeneous  ground  compare  favorably 
with  those  of  previous  theoretical  studies. 

The  more  difficult  cases,  involving  twtvlaycr  earth  models,  generally  require 
less  than  165  sec  of  computer  run  time  for  layer  thicknesses  up  to  10  m.  rod 
lengths  up  to  30  m.  and  wire  lengths  up  to  1(H)  in.  For  these  cases,  best  accur- 
acy is  obtained  with  electrode  radii  less  than  0.01  m.  Specific  problems  such  as 
the  coincidence  of  the  two-layer  interface  with  the  electrode  center  and  subsur- 
face emplacement  of  a horizontal  wue  are  also  discussed 

Tire  computer  program  developed,  written  in  BASIC  computer  language  is 
printed  out  in  full  in  the  appendix,  and  ditticult  cases  involving  lengths  run 
times  and  numerical  inaccuracies  are  listed  at  the  end  ol  the  repoit  No  attempt 
was  made  to  catalog  results  for  variations  in  all  the  pat.  melt  is  that  the  program 
can  consider.  Instead,  it  is  believed  that  the  results  mescnicd  aie  suMicietit  evi- 
dence for  the  capabilities  ol  the  progiant 


IV 


A COMPUTER  PROGRAM  TO  DETERMINE  THE 
RESISTANCE  OF  LONG  WIRES  AND  RODS  TO 
NONHOMOGENEOUS  GROUND 

Steven  A.  Arcone 


INTRODUCTION 

Good  ground  contacts  are  often  necessary  for  the 
protection  of  equipment  and  personnel  against  exces- 
sive electrical  transients  or  overloading.  The  idea  is 
that  when  an  excess  of  current  is  sent  through  the  cir- 
cuit a fault  system  activates  to  divert  this  current 
into  the  earth  through  suitably  placed,  low-resistance 
ground  connections.  The  connections  can  be  a simple 
arrangement  of  a metal  rod  placed  in  the  ground  or  a 
more  complicated  arrangement  of  rods  or  buried 
wires.  In  every  case,  however,  the  resistivity  of  the 
earth  materials  encountered  is  an  important  design 
parameter 

In  the  Arctic,  where  permafrost  or  seasonally  fro- 
zen ground  is  encountered,  a knowledge  of  earth  re- 
sistivity becomes  all  the  more  important  for  engineer- 
ing purposes.  Hoekstra  et  al.  (1975)  and  Scllmann  et 
al.  (1974),  using  electromagnetic,  noncontact  methods 
of  resistivity  surveying,  showed  that  resistivity  in  the 
Arctic  is  highly  variable  but  that,  with  proper  survey- 
ing techniques,  suitable  grounding  sites  can  be  located. 
Their  work  also  showed  tiiat,  in  the  Fairbanks  area, 
values  of  earth  resistivity  are  most  commonly  found 
between  100  and  a few  thousand  ohm-meters  over 
thawed  and  frozen  sediments  of  varying  ice  content 
and  that  this  range  is  oft<m  spanned  in  the  active  layer 
in  the  course  of  a year  at  any  location.  This  is  due  to 
seasonal  thaw  and  changes  in  the  active  layer  depth, 
which  must  also  be  considered  for  grounding  applica- 
tion. 

Previously.  Sunde  (1949)and  Tagg  (1964)  theo- 
retically considered  these  grounding  problems  in  non- 
homogcncous  earth.  However,  Sunde  did  not  pursue 
the  effects  of  variations  in  layer  thickness  or  the 


effects  of  changing  the  penetration  depth  of  a vertical 
rod.  Tagg’s  analysis  is  unsuitable  for  rapid  calculations 
because  he  presented  an  additional  resistance  due  to 
the  penetration  of  the  rou  into  the  lower  earth  layer 
that  had  to  be  added  to  the  resistance  of  the  portion 
of  the  rod  in  the  upper  layer.  After  these  values  were 
found,  he  presented  a penetration  factor  which  had  to 
be  multiplied  by  the  sum  of  these  resistances,  making 
the  calculation  of  resistance  to  ground  laborious. 

Since  calculations  of  this  type  are  so  important  to 
electrical  engineering  in  the  Arctic,  a computer  pro- 
gram was  developed  in  this  study  for  rapidly  calculat- 
ing the  resistance  to  ground  of  either  a vertical  rod  or 
a horizontal  wire,  both  of  specified  length  and  cross 
section,  for  a one-  or  two-layered  earth  model.  The 
program's  accuracy  was  verified  by  comparing  the  re- 
sults with  those  of  Sunde  (1949)  for  the  simpler  cases. 
No  attempt  was  made  to  catalogue  a series  of  earth- 
resistance  curves  for  variations  in  all  available  para- 
meters. The  total  computer  program  is  presented  in 
the  appendix. 


THEORY  OF  EARTH  GROUNDING  FOR 
SIMPLE  ELECTRODES 

Figure  1 shows  an  idealized  case  of  a simple  ground- 
ing configuration.  A hemispherical  electrode  of  radius 
r0  inserted  at  the  surface  of  an  earth  of  resistivity  p. 
measured  in  ohm-meters,  delivers  a current  of  / am- 
peres to  the  ground.  The  homogeneity  of  the  ground 
permits  the  current  to  flow  symmetrically  away  from 
the  electrode  in  hemispherical  fashion  such  tiiat  the 
current  density  J,  measured  in  coulombs/(sec-mJ ). 
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Figure  l.  Hemispherical  electrode  discharging  current  into  a 
homogeneous  earth. 


J* 


Figure  2.  Cylindrical  coordinate  system  on 
a resistively  layered  earth  used  for  solving 
Laplace 's  equation. 


J = 


A 


(1) 


where  r is  the  radial  distance  from  the  electrode  and  ar 
is  the  radial  unit  vector. 

The  electric  field  E,  measured  in  volts/meter,  which 
determines  the  direction  of  current  flow,  is  found 
from  the  relation 


In  cases  where  the  earth  may  be  idealized  as  two  or 
more  uniform  layers  of  different  resistivities,  and  the 
electrode  geometry  becomes  more  complicated,  theo- 
retical solutions  can  only  be  presented  in  integral  or 
series  form  for  which  a computer  must  be  used  to  ob- 
tain numerical  answers. 

The  common  procedure  used  was  to  solve  Laplace’s 
equation  for  the  electrostatic  potential 


-*  -* 
E = pj 


(2)  V2K=0 


(7) 


and  the  electrostatic  potential  V,  measured  in  volts  at 
any  distance  r.,  is  defined  as 

r | 

V=  J E • d?  (3) 

ro 

where  dr  is  the  incremental  radial  distance  of  integra- 
tion and  the  dot  product  is  understood  to  be  taken. 
Carrying  out  the  integration,  we  find  that 


For  unit  current  delivered  to  the  electrode,  the  re- 
sistance R to  uniform  ground  out  to  a distance  rt  is 
then 

/?=-=-£-  (— ohms . (5) 

/ 2n  \r0  rj 

If  the  return  electrode  at  r,  is  considered  at  infinity, 
then 

R = S~  (6) 

2nr0 

is  the  resistance  to  ground  of  the  hemispherical  electrode. 


in  cylindrical  coordinates,  as  illustrated  in  Figure  2.  In 
the  figure,  the  origin  is  taken  at  the  surface  of  the 
earth,  z is  the  depth  below  the  surface,  r is  now  the 
cylindrical  radial  distance,  and  0 is  the  angular  coordi- 
nate. Laplace’s  equation  must  then  be  solved  subject 
to  certain  constraints  upon  the  current  flow  and  the 
potential  itself.  The  constraints  are  that,  at  the  inter- 
face of  any  two  layers,  including  that  between  air  and 
earth,  there  must  be  a continuity  of  potential  V and 
of  normal  current  density  J?.  Since,  within  each  medi- 
um. the  current  density  is  defined  as 

J = Elp  (8) 

where 

E = -VV  (9) 

all  constraints  can  be  mathematically  expressed  in 
terms  of  V.  Air  is  considered  infinitely  resistive,  so 
that  J7  = 0 at  z = 0. 

MATHEMATICAL  PROCEDURE 

The  approach  used  is  based  on  the  division  of 
either  electrode,  the  horizontal  wire  or  vertical 


rod.*  into  discrete  segments.  Considering  one  seg- 
ment at  a time,  its  electrostatic  potential  at  every 
other  segment  position  is  computed.  The  potential 
at  each  segment  is  then  the  sum  of  the  potentials  de- 
rived from  all  the  other  segments.  With  unit  current 
delivered  to  the  entire  rod,  the  resistance  to  ground 
is  then  the  sum  of  all  the  segment  potentials.  This 
process  demands  n 2 calculations  for  n segments.  To 
reduce  this  number,  the  total  potential  is  well  ap- 
proximated by  summing  the  potentials,  calculated  at 
the  electrode  centers,  of  all  other  segments. 

Since  the  potential  developed  around  each  elec- 
trode segment  is  independent  of  the  angular  coordi- 
nate <p  when  the  origin  is  placed  in  the  center  of 
that  segment,  eq  7 becomes 


i±(,rn+*iym0 

r dr  \ dr/  d z2 


(10) 


when  expressed  in  differential  form.  The  solutions 
are  stated  without  proof,  as  the  general  procedures 
for  determining  them  are  given  in  many  texts  (e.g., 
Stratton  1941,  Ward  1967).  All  the  solutions  eval- 
uated at  the  electrode  center  are  of  the  integral 
form 

oo 

V-  J /(X,z, w,z, ,p,  2» A:, )70(Xr) JX  (11) 
o 

where  X = variable  of  integration 

2 = depth  of  the  center  of  the  rod 
n = depth  to  the  center  of  the  segment  in 
question 

2,  = depth  to  the  two-layer  interface 
p,  2 = resistivity  of  either  layer 
= (PrPi/Pi+Pi) 

J0  = zero  order  Bessel  function  of  the  first 
kind 

r = electrode  radius 
cA  = incremental  change  in  X. 

For  the  following  cases,  the  term  segment  refers 
to  that  particular  segment  of  the  electrode  whose 
potential  we  arc  calculating  at  the  electrode  center. 
Four  different  cases  arise,  depending  on  whether 
the  center  of  the  electrode  or  the  center  of  the  seg- 
ment is  in  layer  I or  layer  2. 


* The  distinction  between  rod  and  wire  is  only  on 

the  basis  of  radii. 


Case  1:  Segment  in  layer  1,  electrode  center  in 
layer  1 : 

oo  j 

V = J (a  e^+fie^'+y-  exp(-X|2-n|j 


JJ\r)d\ 


where 


A = - 


p,/  jX,  exp  [-X(22,-n)]-e~x”  ( 
4 n j X,  exp(-2X2,)+  1 ( 

B=  p,/r*,  exp(-2X2,Xe'to+c*"n 
4rr  L X,  cxp(-2Xz  , )+ 1 J 


(12) 

(13) 

(14) 


Case  2:  Segment  in  layer  1.  electrode  center  in 
layer  2: 


V = J Ce'^J0(\r)d\ 


(15) 


where 


p ,/  r An  +/>Art  1 

c = — J-  (i-  x. ) I £ — — . (16) 

47t  |_X,  exp(-2Xz,)+ 1J 

Case  3:  Segment  in  layer  2,  electrode  center  in 
layer  1 : 

oo 

V = J Die'**  +e^2)yo(Xr)JX  (17) 


where 


4 Tt 


*,)[ — 1.  (18) 

|_X , exp(-2X2,)+  1 J 


Case  4:  Segment  in  layer  2,  electrode  center  in 
layer  2: 

oo 

V = J H exp(-X|2-«|  ) + /■.'  e'^Jy^XTyX 
o 

(19) 


where 


4tt  ( [kx  exp(-2XZj )+l J 

- cxp(2Xr,  )|. 


(20) 


Resistance  to  ground  as  a function  of  electrode  length  for  Figure  4.  Resistance  to  ground  as  a function  of  electrode  length 

•od  placed  in  various  one - and  two-layer  earth  models.  horizontal  wire  lying  on  various  one-  and  two- layer  earth  models. 


0 2 4 6 8 10 


Z,  (m) 

Figure  5.  Resistance  to  ground  of  a vertically  driven 
rod  as  a function  of  first-layer  thickness  z. 

Tliese  equations  govern  both  horizontal  and  verti- 
cal electrode  segments.  For  a horizontal  electrode, 
either  at  the  surface  or  buried,  z = n.  When  a verti- 
cal electrode  penetrates  both  layers,  the  following 
assumptions  are  used: 

1.  The  current  per  unit  length  V,  in  the  portion  of 
the  electrode  in  the  first  layer,  and  the  current  per  unit 
length  Kj  in  the  portion  of  the  electrode  in  the  second 
layer,  are  related  to  the  total  current  / by  the  equation 
(Tagg  1064) 

Vi  + yV'2  = / (21) 

where  l.  , and  f2  are  the  electrode  lengths  in  the 
first  and  second  layers,  respectively. 

2.  T,  and  Y2  arc  related  by  the  equation  (Tagg 
1064) 


The  computer  program  developed,  called  RHS1ST. 
is  written  in  BASIC  computer  language,  and  is  listed 
in  the  appendix.  The  program  performs  the  segmen- 
tation of  the  electrodes,  integrates  eq  12,  15.  17  and 
10,  and  sums  the  resulting  potentials  to  find  tire  total 
resistance  to  ground  for  unit  applied  current.  Hither 
homogeneous  ground  or  two-layer  ground  models  can 
be  considered.  A definition  of  all  the  parameters 
used  is  listed  at  the  start  of  the  program. 

RESULTS 

To  reduce  the  number  of  cases,  only  the  following 
values  were  used:  resistivity  values  of  100  and  1000 
ohm-m.  vertical  rod  radius  of  0.05  m.  and  horizontal 
wire  radius  of  0.00125  m.  Tliese  electrode  dimensions 
were  chosen  to  enable  a comparison  with  Sunde’s 
(1949)  results  for  homogeneous  ground  to  be  made. 

Tire  program  considers  the  electrodes  to  be  only  at  or 
below  the  surface. 

The  resistances  to  ground  of  the  two  types  of 
electrodes,  as  a function  of  increasing  electrode  length 
for  various  one-  and  two-layer  cases,  arc  plotted  in 
Figures  3 and  4.  The  dashed  curves  are  Sunde’s  re- 
sults and  should  be  compared  with  the  100-ohm-m 
solid  curves  that  are  nearest  them.  Tliese  favorable 
comparisons  ensure  the  validity  of  the  program  for 
wires  less  than  100  m long.  The  two-laver  curves 
show  the  convergence  of  the  resistance  to  ground 
of  the  one-layer  and  two-layer  models  for  long  rod 
lengths.  The  curves  reveal  that  the  horizontal  wire 
gives  superior  grounding  per  meter  of  length  when 
the  upper  layer  is  more  conductive  (100  ohm-m/ 1 000 
ohm-m  case)  than  the  lower  layer,  whereas  the  driven 
rod  is  superior  when  the  lower  layer  is  more  conduc- 
tive (1000  ohm-m/100  ohm-m  case)  than  the  upper 
layer.  This  is  logical,  since  the  better  performance 
is  exhibited  by  the  electrode  maintaining  more  con- 
tact with  the  100-ohm-m  earth. 

Tlie  variations  in  resistance  to  ground  as  a func- 
tion of  first-layer  thickness  for  the  two  different 
electrodes  are  plotted  in  Figures  5 and  6.  This  is 
an  important  case  to  consider  for  frozen  ground  ap- 
plications where  the  active  layer  varies  in  thickness 
during  the  year.  In  areas  of  seasonal  frost,  where 
only  the  upper  layer  experiences  freezing,  the  30-m- 
long  vertical  rod  electrode  is  superior  to  the  100-m- 
long  wire  for  all  layer  thicknesses,  as  can  be  seen  by 


5 
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Figure  6.  Resistance  to  ground  of  a horizontal  wire 
resting  on  the  surface  as  a function  of  first-layer 
thickness  z. 

comparing  the  broken  curves  of  the  two  figures.  In 
permafrost  regions,  the  100-m-long  wire  is  superior 
when  the  active  layer  is  thawed,  as  can  be  seen  by 
comparing  the  solid  curves  of  the  two  figures. 

These  results  vary  with  electrode  length,  but  they 
indicate  the  value  of  the  program  for  considering 
these  important  situations. 


division,  which  works  to  a tolerance  of  about  5% 
until  the  radius  equals  0.01  m;  radii  smaller  than 
this  may  therefore  require  more  segments.  Thus, 
it  is  recommended  that  vertical  electrode  radii  be 
greater  than  I cm  to  save  on  computer  editing  and 
run  time. 

2.  Center  of  the  vertical  electrode  (rod)  near 
the  two-layer  interface.  This  method  of  segmenting 
the  electrode  requires  that  each  segment  lie  entirely 
within  one  layer.  Therefore,  if  a desired  electrode 
is  found  to  have  a center  coincident  with  the 
modeled  interface,  it  is  recommended  that  it  be 
lengthened  a distance  equal  to  six  radii  or  lowered 

a distance  of  three  radii  to  avoid  this  situation. 

3.  Horizontal  electrodes  (wires)  buried  beneath 
the  surface.  To  check  on  the  accuracy  of  the  results 
for  these  cases,  it  is  recommended  that  the  electro- 
static potential  be  recalculated  for  all  segments  where 
r > 1 m by  decreasing  d\  by  a factor  of  10  or 
greater. 

CONCLUSIONS 

The  program  RESIST  is  an  effective  means  of  cal- 
culating the  resistance  to  ground  of  simple  electrodes. 
Tlie  number  of  variables  allows  a great  deal  of  flexi- 
bility in  dealing  with  two-layer  earth  models.  Use  of 
this  program  is  facilitated  by  interactive  aids  whereby 
the  user  inputs  data  in  response  to  formatted  ques- 
tions. When  specific  applications  for  arctic  regions 
are  concerned,  seasonal  resistivity  and  depth  changes 
must  be  obtained  for  the  site  selected  before  a spe- 
cific electrode  type  is  decided  upon.  Then,  curves 
such  as  those  in  Figures  5 and  6 can  be  used  to 
determine  the  seasonal  change  of  the  ground  resis- 
tance for  specific  electrode  types  so  that  a year-long 
grounding  system  can  be  designed. 


PROBLEM  AREAS 

The  following  situations  should  be  avoided  when 
using  the  computer  program  listed  in  the  appendix: 

I.  Vertical  electrode  (rod)  radii  less  than  0.01  m. 
The  program  approximates  each  segment  of  the 
electrode  by  its  midpoint.  This  may  become  inac- 
curate with  the  present  algorithm  for  electrode 
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APPENDIX:  COMPUTER  PROGRAM 
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RESIST  (continued ) 


0430  ' AND  'S-2.MED2*  ONLY)  - INDICATES  TU  THE  USER  AND  TO  THE  COM-  10930  LET  T=T- . 00054* ( 3/Y ) ~4- . 00029333* < 3/Y > "3+ . 0001 3558* ( 3/Y ) 

0440  ' PUTER  WHICH  SHORTCUTTING  PATH  IS  FOLLOWED  WITHIN  A SUBPROGRAM  10940  LET  J=Y~ ( - . 5 ) *F*COS ( T ) 

,0430  ' (THESE  SHORTCUTS  ARE  DESIGNED  TO  SAVE  COMPUTER  TIME  AND  SPEED  10950  SUBEND 


<30  ' SPECIAL  CASE  WHERE  DOTH  THE  CURRENT  SOURCE  AND  THE  REFERENCE  11930  LET  C»  1 / ( K1  »EXP(  -2*L«Z1  ) 1 1 ) 

<40  ’ POINT  ARE  AT  THE  SURFACE.  THE  EOUATION  IS  REDUCED  TO  A SERIES  11940  LET  K=K9»C*FNH ( Z > * J 

450  • EXPANSION  LIKE  THAT  ON  PAGE  51  (EO.  2.39)  IN  SUNDE  (1949).  11950  IF  APS < ENH( Z ) ) IE-3  THEN  12000 
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